Thermodynamics of carrier-mediated magnetism in semiconductors 
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We propose a model of carrier-mediated ferromagnetism in semiconductors that accounts for the 
temperature dependence of the carriers. The model permits analysis of the thermodynamic stability 
of competing magnetic states, opening the door to the construction of magnetic phase diagrams. 
As an example we analyze the stability of a possible reentrant ferromagnetic semiconductor, in 
which increasing temperature leads to an increased carrier density, such that the enhanced exchange 
coupling between magnetic impurities results in the onset of ferromagnetism as temperature is raised. 



Ferromagnetic semiconductors (FS) show important 
and potentially useful differences from their metallic 
counterparts. For example, if the magnetism in a mag- 
netically doped semiconductor is mediated by carriers, 
then changes in the carrier density induced by light or 
applied bias may significantly alter the the exchange in- 
teraction between the carriers and magnetic impurities. 
When this effect is sufficient to turn the ferromagnetism 
on or off, there arise intriguing possibilities for light- or 
bias-controlled ferromagnetism [l], 0, [H, 0] not possible in 
conventional ferromagnetic metals. 

Another approach to controlling ferromagnetism in FS 
materials is to exploit the strong temperature depen- 
dence of carrier density that is the hallmark of semi- 
conductors. Despite this dependence, most theoretical 
analysis of FS has assumed a metallic-like picture [B|, Q 
in which carrier density is treated as independent of tem- 
perature. Consequently, a number of experimental fea- 
tures observed in certain FS materials remain incom- 
pletely explained. Examples include the metal-insulator 
transition [7| in Mn-doped GaAs: the role of the impu- 
rity band in oxide FS materials [a, 0] ; and the resistivity 
peak in Mn-doped GaAs observed near the Curie tem- 
perature, T(7, which is usually attributed to temperature- 
dependent scattering [l(| just as in metals 11 1. In each of 
these cases, correct treatment of the temperature depen- 
dence of the carrier density could substantively change 
our understanding of the observed phenomena. 

In this paper we develop a theoretical model of ferro- 
magnetism in semiconductors that includes the temper- 
ature dependence of the carriers. By providing a way to 
analyze the thermodynamic stability of different compet- 
ing magnetic states, this model makes possible the self- 
consistent calculation of the magnetic phase diagram of a 
FS. Here we use the model to calculate the temperature- 
dependent magnetization of a simple generic FS. In con- 
trast to the standard monotonic decay of the magneti- 
zation with increasing temperature found in metals, we 
demonstrate the possibility of stable "reentrant" ferro- 
magnetism in semiconductors: as the temperature is in- 
creased the resulting higher density of thermally excited 
carriers can enhance the exchange coupling between mag- 



netic impurities — and thereby increase the magnetization 
over some range of temperatures. Of course, whether 
such a possibility can be realized even in principle de- 
pends on the thermodynamic stability of the reentrant 
phase relative to other magnetic states. To properly an- 
alyze this competition requires a theoretical framework 
for computing the free energy of each possible magnetic 
state. Our model provides this framework. 

We begin by considering a FS doped with donors of 
density Nd and with magnetic impurities of density iVj. 
For simplicity we assume that no acceptors are present, 
that electrons are the only carriers, and that the mag- 
netic impurities are electrically neutral; the more general 
case is straightforward. The interaction between electron 
spins Si and localized impurity spins Jj is taken to be 



ex/]S(ri - Rj)si 
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where T ex is the exchange coupling and ?i (Rj) is the 
position of the carrier (impurity) . We assume a nondc- 
generate semiconductor in which the conduction band is 
separated from the donor level (or impurity band) by Ed 
in the absence of magnetic order, as in Fig. [Ha). When 
magnetic order is present the conduction band and donor 
level experience spin splittings of A [l^] and A^ [l]| , re- 
spectively, as in Fig.QJb). For simplicity we assume that 
the ratio A^/A = 7 has a fixed, material-specific value. 

These spin splittings can arise from either an applied 
magnetic field or from the carrier-impurity interaction, 
and in general will modify the temperature-dependence 
of the electron density [14| ■ We can obtain a expression 
for the density n of conduction electrons as follows. Elec- 
troneutrality requires that n + N^ = Nd, where N® is the 
density of neutral (non-ionized) donors. If we take the 
effective density of states in the conduction band as N c = 
{\/A)(2m*k B T/TTh 2 f/ 2 d5], then the electroneutrality 
condition can be satisfied by introducing the chemical po- 
tential /i that satisfies N c exp(^i/fcsT) cosh(A/2fcsT) = 
N d / [1 + 2 exp [(n + e d )/k B T] cosh(A d /2fc B T)]. This is 
equivalent to a quadratic equation in the electron den- 
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FIG. 1: Schematic energy-band structure, and temperature 
evolution of free energy, (a) In the absence of magnetization 
(M = 0) both the conduction band and the donor level (at 
—Ed relative to the band edge) are unsplit. (b) The onset 
of magnetic order (M 7^ 0) leads to spin splitting A of the 
conduction band and Ad of the donor level . Thick arrows 
represent two spin projections, (c) Free energy T versus aver- 
age electron spin s. Materials parameters were chosen based 
on Gd-doped EuO [Hj]: T ex = 40 eV-A 3 , \e d \ = 20 meV, 
N d = 3.5 x 10 19 cm -3 , m* = 2, lattice constant a = 5.15 A, 
Ni = 4/a 3 ,, 7 = 1, a B = 8 A, and J = 7/2. Behavior at 40 
and 150 K reveals ferromagnetic states. 



sity, and has the solution 
1 



n(s,T)= -N dK (s,T) y/l + 4/K(s,T)-l 



(2) 



Here we have expressed the dependence of n on the split- 
ting A equivalently as a dependence on the average spin 
of the conduction electrons, s = (1/2) tanh(A/2fcsT). 
We have also defined the auxiliary function k(s,T) = 
(N c /N d ) exp(- £d /fc B T)(l - 4 S 2 )^" 1 )/2/[(l + 2s)^ + (1 - 
2s) 7 ]. In addition, for later use it is here convenient 
to introduce the relative density of conduction electrons 
v(s,T) = n{s,T)/N d , so that < u{s,T) < 1. 

In the magnetically ordered state, both the conduction 
electrons and the electrons bound to donors experience 
the field from the ordered impurity spins. The impurity 
spins in turn experience both a long-ranged uniform field 
from the delocalized conduction electrons and a short- 
ranged field from the bound donor electrons. This situ- 
ation can be described using two coupled order param- 
eters: the average spin s of the conduction electrons as 
defined above; and the average spin m = M/Nigifis of 
the impurities (here M is the total magnetization due to 
the impurities and gt is the electron ^-factor). 

Using these two order parameters, we can write the 
free-energy density of the system in the form 

T = Ti + T e - T ex NiN d [v(s, T)s + v d (s, T)s d ] to. (3) 

Here Ti = Ti (to) is the contribution from the entropy of 
the impurity spins, and T e = T e (s) is the contribution 
from the entropy of the conduction electrons. The for- 
mer can be obtained by expressing the free energy of the 



impurity spins as a function of an external magnetic field 
and then applying a Legendre transformation [161 ]. This 
gives 

sinh[(2J+ l)a(m)]" 



2Ja(m) — In ■ 



sinh [a (to)] 

(4) 

1 (x) denotes the 



where a(m) = B~j (m/J)/2J, and Bj 
inverse of the Brillouin function (lBj . 

The electron contribution J- e (s) can be derived using 
a similar approach. We first define the density- weighted 
average spin of conduction and donor electrons, 



t(s, T) = v(s, T)s + [1 - v(s, T)]s d , 



(5) 



where the average spin of donor electrons is given by 
s d = (1/2)[(1 + 2sY> - (1 - 2s) 7 ]/[(l + 2s)t + (1 - 2s)^. 
We next invert the function r(s, T) by solving Eq. (O for 
s = s(t, T). This function can then be used to obtain 
the free energy of the conduction and donor electrons: 

Ms) = k B TN d J o ln T ±^j^. (6) 

The third term in Eq. ([3]) is the mean-field approxi- 
mation for the internal energy described by the exchange 
Hamiltonian of Eq. (JTJ) . This term represents the cou- 
pling of the order parameters s and to. The expression in 
square brackets makes explicit the separate contributions 
from conduction electrons and donor electrons. Since 
the conduction electrons are delocalized, they mediate a 
long-ranged interaction between impurity spins. In con- 
trast, the interaction between donor electrons and impu- 
rity spins is short-ranged, and is controlled by a contact 
term given by the value of the donor wavefunction at the 
impurity site. We turn now to evaluating this interaction 
by deriving an explicit expression for the term v d (s,T) 
appearing in Eq. ([3]). 

We use a "two-color" percolation model to represent 
the short-ranged interaction between the randomly dis- 
tributed donor electrons and impurity spins (ItI [l8j . This 
model is a generalization of a more standard percolation 
approach originally proposed for dilute ferromagnets Il9( 
and recently applied to magnetic semiconductors [20(. 
Within the two-color model, an interaction is counted 
for each pair of sites whose spatial separation is less than 
R c = [S c /(47r/3)] 1 / 3 (iV l 7V0)- 1 /^ where B c w 2.7 is the 
average coordination number [171 ]. The density of such 
pairs within the infinite percolation network is N c iB c /2, 
where N d = 2/(4ti\R c /3) 3 is the average density of donor- 
and impurity-sites belonging to the network. 

Within the mean-field approximation, the contribu- 
tion of these pairs to the internal energy can be shown 
to be Ud = -r e2 |V(0)| 2 exp(-2i? c / ai3 ) {l/2)N cl B c ms dl 
where "0(0) is the value of the donor wavefunction at the 
origin (|t/;(0)| 2 = l/ira 3 B for hydrogenic donors). Com- 
paring this result to Eq. ([3]) implies that 

v d = |0(O)| 2 (iV l 7V (i )- 1 / 2 (l-^) 1/2 exp(-2i? c /a B ). (7) 
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FIG. 2: Temperature dependence of order parameters s and 
to, showing reentrant behavior at T — Tci as temperature is 
raised. Curie temperatures Tbi,2,3, at which to and s both 
vanish, define phase boundaries between two ferromagnetic 
and two paramagnetic states, (a) Inset: Temperature depen- 
dence of the combined average spin from all electrons. Mate- 
rials parameters are the same as in Fig. [JJ 



Having now obtained explicit expressions for all the 
terms in Eq. ([3]), we can directly obtain the temperature 
dependence of the order parameters s and m by minimiz- 
ing the free energy T . This results in 



1 



■ tanh 



T ex N i md[i){s 1 T)s] 



ni = JBj 



2k B T dr{s,T) 



JT ex N d v(s,T)s 
k B T 



(8) 



(9) 



where v(s,T) = v(s,T) + Vd(s,T)sd/ s. By expanding 
Eqs. ([8]) and ([9]) for small s and m we obtain an implicit 
expression for the critical temperature: 

T c = T c *(0,T C )/ h + (1 - 7>(0, T c )f 2 , (10) 

where T c = {T ex /k B )[NiN d J{J + 1)/12] 1 / 2 is the Curie 
temperature in the limit of completely ionized donors Q . 

We now turn to exploring the predictions of our model 
for a realistic example. We choose materials parameters 
approximately corresponding to Gd-doped EuO [2l|, a 
magnetic material known to exhibit strong temperature 
dependence of the carrier density. Fig. [3c) shows the re- 
sulting free energy J 7 as a function of conduction-electron 
spin s, at the temperatures 40, 60, 150, and 300 K. A fer- 
romagnetic state (at 150 K) appears at higher tempera- 
ture than a paramagnetic state (at 60 K), and is thus 
a reentrant ferromagnetic state. This phenomenon is a 
direct consequence of the increased number of thermally 
excited carriers, which sufficiently increase the exchange 
coupling between magnetic impurities to overcome the 
additional entropic cost of a magnetically ordered state. 



Figure [2] shows the temperature dependence of s(T) 
and m(T) for the same materials parameters. For this 
example there are three different solutions of Eq. (jTUJ) 
for Tc, as shown. The reentrant ferromagnetic state ex- 
ists in the temperature range Tci < T < Tci- The in- 
set shows the temperature dependence of the combined 
average spin from all electrons (conduction, donor, and 
impurity) , revealing behavior very different from the con- 
ventional monotonic decay [l5j |. 

In the standard theoretical description of dilute mag- 
netic semiconductors the ferromagnetism is mediated en- 
tirely by "free" carriers (usually holes in the valence 
band) 0, Q . An important check of our theory is that it 
reproduces earlier results obtained in this limit. We thus 
consider the case with all donors ionized, n = Nd- In this 
regime the entropy contribution T e (s) must be replaced 
with its degenerate counterpart for the internal energy, 
JF dc s( s ) = (3 n £ F /l0) [(1 + 2s) 5 / 3 + (1 - 2s) 5 / 3 ] , where 
E F = (n 2 /2m*)(37r 2 n) 2 / 3 is the Fermi energy. Sub- 
stituting this expression into Eq. ([3]) and expanding 
T in s and m near the Curie temperature, we obtain 
T c = {T 2 ex /k B )NiN J{J + 1)/12, where N = 3n/2E F is 
the electronic density of states at Ep . This is the stan- 
dard expression for Tc given, for example, as Eq. (7) of 
Ref. U. 

The possibility of reentrant ferromagnetism in semi- 
conductors was first discussed over forty years ago [22l |. 
and again recently [23] . In neither case was the thermo- 
dynamic stability of the magnetism discussed. Moreover, 
the role of the donor electrons was not included, an omis- 
sion that can lead to contradictions. For example, in Ref. 
[23| it was assumed that the conduction-electron density 
has no explicit dependence on the magnetization. This 
assumption is justified only if the spin-splitting of the 
donor level (Ad) and conduction band (A) are equal, i.e. 
7 = 1. However, the self-consistency equation given in 
Ref. [23l | for the magnetization of the impurity spins is 
only correct for the case 7 = 0. This can be seen by 
substituting our Eq. (|8|) into Eq. ((9]) and comparing to 
the corresponding equation in Ref. [23]. This internal 
inconsistency has substantial consequences. Specifically, 
for the case 7=1, expanding the self-consistency equa- 
tion for small magnetization in the vicinity of the critical 
temperature yields Tc = T c u(0,Tc). But for the case 
7 = 0, a similar expansion yields the qualitatively differ- 
ent dependence T c = T cy /v{Q,T c ). 

In its present form, our model does not include several 
physical effects which could alter the magnetic ordering 
at low temperature. (1) Hopping and Coulomb fluctua- 
tions will broaden the donor levels into an impurity band 
of localized states [24j]. (2) The neutral donors may form 
bound magnetic polarons (25| by aligning nearby impu- 
rity spins, thereby increasing the donor ionization energy 
Ed- (3) There could be a substantial direct antiferro- 
magnetic interaction between nearest-neighbor impurity 
spins. This would have two consequences: J- in Eq. ([3]) 
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FIG. 3: Temperature dependence of order parameters s and m 
when percolation is omitted (see text). Materials parameters 
are the same as in Fig. [TJ except for m*—2.5. 



would acquire a contribution oc to 2 , while the argument 
in the mean-held expression for m in Eq. (|9|) would be 
reduced by a term oc to. 

All three of these effects would suppress ferromag- 
netism at low temperatures. While the situation is fairly 
complex, we can describe it qualitatively by omitting the 
percolation term in Eq. ([3]), i.e. by setting Vd(s,T) = 
or, equivalently, by taking as — > in Eq. ([7]). Figure 
[3] shows the resulting behavior for s(T) and m(T). The 
striking feature is that the ferromagnetism is absent at 
both low and high temperatures, corresponding to anti- 
ferromagnetic and paramagnetic states, respectively. A 
similar interpretation of reentrant ferromagnetism was 
proposed in Ref. 22( to explain experimental data on 
Furthermore, both of the scenarios for 



(Li,Mn)Se [26 

reentrant ferromagnetism shown in Figs. [Hand [3] are con- 
sistent with the recent experiments in (In,Mn)Se [27j |. 
With the change in Mn-concentration there is a change 
in the number of peaks in the temperature dependence of 
dynamic magnetic susceptibility supporting the existence 
of either two or three distinct critical temperatures (re- 
call Figs. [2] and [3j. While, for simplicity, we have focused 
on the parameters for Eu-based semiconductors, there is 
a need to explore other FS for a possible reentrant ferro- 
magnetism. 

Predictions of reentrant ferromagnetism could be di- 
rectly tested in transport experiments. For exam- 
ple, in conventional (non-reentrant) spin light-emitting 
diodes [4J the emitted light is circularly polarized if the 
injected carriers are spin-polarized. As the temperature 
is raised, there is a monotonic decrease in the degree of 
polarization of the electroluminescence, due to the de- 
creased polarization of the carriers. By using a reentrant 
FS as the spin injector, a nonmonotonic temperature de- 



pendence of the electroluminescence, similar to that for 
s{T) in Figs. [2] or [3j would be expected. Reentrant be- 
havior could also be detected electrically, without a spin- 
LED, by demonstrating non- monotonic temperature de- 
pendence of the magnetoresistance in a semiconductor 
heterojunction 1J, |28| that includes a reentrant FS. 
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